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Abstract: Cephalopod taxonomy is still uncertain, and little is known of the phylogeny of Recent taxa. Biochemical and molecular characters are com¬ 
plementary to morphology, and allow an additional insight into the phylogenetic relationships among cephalopods. Eye lens protein electrophoresis and 
immunological approaches yield data in agreement with traditional taxonomic grouping, but are less suitable for establishing phylogenetic relationships. 
Molecular tools, e. g. the 3’ end of the 16S rDNA gene, have failed to resolve the phylogeny at the suprafamilial level, but seem appropriate at lower levels. 
DNA sequence comparisons (% substitution) show that a direct relationship between taxonomic rank and nucleotide divergence cannot be established, as the 
nucleotide divergence level differs from one taxa to the other. 

Electrophoretic and immunological analyses of eye lens proteins as well as molecular results suggest that sepiolids should be separated from 

Sepioidea. 


Phylogenetic analyses allow hypotheses of biologi¬ 
cal evolution based on various criteria to be tested. 
Morphological classification of taxa has, historically, been 
the first approach considered, and the analyses of anatomi¬ 
cal structures are still of prime importance in the construc¬ 
tion of hypotheses on the evolutionary history of organ¬ 
isms. But more recently, the development of biochemical 
and molecular techniques, that allow an insight into genetic 
structure, have opened new perspectives in taxonomy and 
phylogeny. Paleontology provides, however, the only direct 
means of calibrating evolutionary events. Both paleontol¬ 
ogy and embryology allow recognition of homologous 
characters derived from a common ancestor, the ancestors 
being deduced from the characters of terminal groups. 

As far as cephalopods are concerned, the taxonomic 
status of many groups is still uncertain, and little is known 
of the phylogeny of Recent taxa. The main reason for this is 
that although the fossil record contains a wealth of ecto- 
cochleate forms (all extinct except Nautilus ), this is not the 
case for Recent, mainly soft-bodied animals with a reduced 
or absent internal shell. So far only a few authors have con¬ 
sidered the phylogenetic systematics of the group 
Cephalopoda. Fioroni (1981) was one of the first to use the 
Hennig (1966) approach to systematics. Berthold and 
Engeser (1987) established adelphotaxonomic relationships 
by the identification of synapomorphies in fossil and 
Recent taxa and proposed a phylogenetic classification of 
35 subordinate taxa. The cladograms presented in these 
papers are not totally in agreement, drawing attention to the 
necessity of using additional criteria, such as biochemical 


and molecular, for taxonomic and phylogenetic reconstruc¬ 
tion of cephalopod systematics. 

All approaches, morphological, biochemical, and 
molecular, are complementary but each has advantages and 
disadvantages. Molecular analyses provide new tools to test 
the hypotheses based on morphology, and help to reformu¬ 
late them in some cases. Morphological characters are more 
accessible, easier and less costly to analyze, allow compari¬ 
son of extant and fossil forms, but can be more subjective. 
Biochemical and molecular characters are more objective, 
are potentially very abundant, but are not always easy to 
analyze. For instance, the fact that the nucleotides at each 
position can exist in only four states could be an important 
source of homoplasy. It is also well known that different 
portions of the same gene have not all the same probability 
of variation, and the rate of evolution can be different for 
the same gene in different taxa. 

This paper is aimed to compare the results obtained 
using biochemical (electrophoretic and immunological) and 
molecular techniques for phylogenetic reconstruction of 
cephalopods using a number of species, comprising 
octopods and decapods. 

MATERIALS AND METHODS 

Fresh cephalopod tissue samples were obtained 
from various sources (Table 1). They were either stored at 
-20°C prior to electrophoretic and immunological analyses, 
or alcohol preserved for nucleotide sequencing. 
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Table 1. List of cephalopod species analyzed by eye lens protein electrophoresis and mtDNA sequenc¬ 
ing. Their geographical origin is indicated as well as the source of the data: 1, Bonnaud et al, 1994; 2, 
Tranvouez and Boucher-Rodoni, 1990; X, present paper; (-), no data. 


Origin 

Species 

Eye lens 

mtDNA 

East Atlantic (Biscay) 

Sepia officinalis Linne, 1758 

X 

1 

Mediterranean (Banyuls) 

S. orbignyana Ferussac, 1826 

2 

1 

SW Pacific (New Caledonia) 

S. latimanus Quoy and Gaimard, 1832 

X 

1 

English Channel (Roscoff) 

Loligo vulgaris Lamarck, 1798 

2 

1 

SW Pacific (New Caledonia) 

Sepioteuthis lessoniana Lesson, 1830 

X 

1 

English Channel (Roscoff) 

Sepiola atlantica Orbigny, 1840 

X 

1 

Mediterranean (Banyuls) 

Rossia macrosoma (Delle Chiaje, 1829) 

X 

1 

Mediterranean (Banyuls) 

Todaropsis sp. 

X 

(-) 

Pacific Ocean (Hawaii) 

Todarodes sp. 

(-) 

1 

Mediterranean (Banyuls) 

Octopus vulgaris Cuvier, 1797 

2 

(-) 

Mediterranean (Banyuls) 

Eledone cirrhosa (Lamarck, 1798) 

2 

1 

SW Pacific (New Caledonia) 

O. cyanea Gray, 1849 

X 

1 

SW Pacific (New Caledonia) 

O. sp. 

X 

(-) 

East Atlantic (Mauritania) 

Graneledone verrucosa (Verrill, 1881) 

X 

(-) 

East Atlantic (Mauritania) 

Opistoteuthis agassizii (Verrill, 1883) 

X 

(-) 


For biochemical approaches (electrophoresis and 
immunology), the protein chosen should be stable enough 
that individual physiological changes do not influence the 
observed differences between species, but it should be vari¬ 
able enough to reflect taxonomic and phylogenetic differ¬ 
ences. Eye lens proteins and hemocyanin were tested here. 

Eye lens proteins were extracted according to the 
protocols described in Tranvouez and Boucher-Rodoni 
(1990), and analyzed on precast polyacrylamide gels 
(ExcelGel SDS Gradient 8-18%, Pharmacia). A band pres¬ 
ence/absence matrix was computed and processed by the 
NTSYS-pc program (Rohlf, 1990), using Sahn clustering 
(Sneath and Sokal, 1973) with UPGMA and Neighbor- 
Joining (NJ) methods (Saitou and Nei, 1987), to produce 
phenograms. PAUP 3.1 (Swofford, 1985) was used to esti¬ 
mate phylogenetic trees. 

Protein antigenic properties are supposed to allow 
the estimation of immunological distance between taxa 
(Tsusumi et al , 1989). The ELISA immunological tech¬ 
nique, an enzyme-linked immunosorbent assay, was adapt¬ 
ed to cephalopod eye lens protein (Boucher-Rodoni et al ., 
1995). Four eye lens antisera were available ( Sepia offici¬ 
nalis , S. orbignyana, Loligo vulgaris , Octopus vulgaris). A 
homologous standard inhibition curve was determined with 
2500-fold diluted serum, and the affinity of heterologous 
samples was then tested. 

The antigenic properties of hemocyanin, the large 
respiratory protein found in the blood of all cephalopods, 
was also used here to estimate taxonomic relationships. The 
immunological distance of various cephalopod species was 
estimated by heterologous reaction against commercial 
keyhole limpet hemocyanin (KLH) and compared to pre¬ 
liminary results obtained with homologous antiserum. 

For molecular analyses, attention was first focused 


on mitochondrial DNA (mtDNA) because of its diversity, 
and because data on various groups, including non-mollus- 
can invertebrates, are well known. mtDNA is a small circu¬ 
lar DNA molecule present in many copies in the mitochon¬ 
dria, it is maternally inherited, and there is no recombina¬ 
tion. The protein genes and the r-RNA genes have both a 
mosaic structure of conserved and variable regions, which 
should allow phylogenetic relationships at various hierar¬ 
chical levels to be analyzed. Nucleotide sequence data from 
the y end of the 16S rDNA gene have already been used to 
analyze phylogenetic relationships among decapod 
cephalopods (Bonnaud et al ., 1994). To determine whether 
the different taxonomic hierarchical levels can be related to 
a given nucleotide percentage of divergence, some 
sequences from two populations of the same species, from 
different species of the same genus, and from different gen¬ 
era, were compared and analyzed in terms of molecular 
divergence (% substitution). 

RESULTS 

The results of eye lens protein electrophoresis of 14 
cephalopod species analyzed by UPGMA and Neighbor- 
Joining methods show that with both methods the taxonom¬ 
ic grouping is respected (Fig. 1). Octopods are always 
grouped together, but are not distinctly separated from 
decapods. Incirrate octopods are a sister group of cirrate 
octopods, but at lower taxonomic levels the genus Octopus 
is not homogeneous. As far as decapods are concerned, the 
relationships among myopsids, oegopsids, and sepioids are 
not clearly defined by either approach. The sepiolids occu¬ 
py a particular position in the NJ analysis, branching as a 
sister group of all other coleoids. The strict consensus tree 
derived from the phylogenetic analysis (PAUP 3.1) is less 
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I- Octopus vulgaris 
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Fig. 1. Phylogenetic tree obtained by the Neighbor-Joining method on SDS-Page electrophoretic analysis of eye lens proteins. 


resolved for octopods (Fig. 2). 

When using protein eye lens immunological proper¬ 
ties to estimate immunological distances for phylogenetic 
analyses, one of the main problems encountered was that in 
some cases the distances were not symmetrical. Table 2 
shows that when comparing the heterologous reaction 
between two species using each species extract in turn as 
antiserum and antigen, the immunological distance is not 
necessarily the same. 

The hemocyanin is currently used as an immuno¬ 
genic agent. In cephalopods it is composed of seven func¬ 
tional units in octopods and Nautilus , and eight functional 
units in decapods (Van Holde et al ., 1992). Fig. 3 shows the 


immunoreactivity of the hemocyanin of various cephalopod 
species against KLH antiserum, estimated by the ELISA 
technique. Nautilus reactivity was the closest to KLH 
homologous reaction, most of the other species being 
grouped together at a rather distant position, except Sepia 
which displayed a reactivity stronger than all other 
coleoids. Preliminary assays with homologous hemocyanin 
antiserum indicate that the distances are readable at high 
hierarchical levels (/. e. distant taxa), but do not discrimi¬ 
nate species. The distance between Sepia and Nautilus was 
the smallest, but the difference was not as important as with 
KLH. And again, as with eye lens proteins, the results were 
not symmetrical. 


Graneledone verrucosa 
Octopus vulgaris 
O. sp. 

O. cyanea 
Eledone cirrhosa 
Opistoteuthis agassizii 
Todaropsis sp. 

Sepia orbignyana 
S. officinalis 
S. latimanus 
Sepiola atlantica 
Rossia macrosoma 
Loligo vulgaris 
Sepioteuthis lessoniana 


Fig. 2. Strict consensus tree (unrooted) of 100 trees obtained by heuristic search (PAUP 3.1; MULPARS option). 
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Table 2. Results of immunological heterologous reaction (ELISA tech¬ 
nique) between species pairs using each species extract in turn as anti¬ 
serum and antigen. Values represent immunological distances estimated 
by optical density differences. 


ANTISERUM 

Sepia officinalis 

<£ 

4 

ANTIGEN 

S. orbignyana 


4 

S. officinalis 


10 

Loligo vulgaris 


15 

L. vulgaris 


17 

S. orbignyana 




13 

Octopus vulgaris 


16 

S. officinalis 


25 

0. vulgaris 


15 

S. orbignyana 


17 

0. vulgaris 


24 

L. vulgaris 


25 


To analyze the relationships of cephalopods at a 
perispecific level (population, species), many authors have 
studied enzymatic polymorphism which remains a valuable 
tool for genetic population analyses (Levy et al. , 1988; 
Carvalho et al., 1992; Brierley et al., 1993, 1995). 

If mtDNA is considered, preliminary results on six 
different populations of Sepia officinalis indicate that nei¬ 


ther 16S, nor cytochrome oxidase Coll or Colli, displays 
adequate variability for taxonomic purposes, whereas enzy¬ 
matic polymorphism does (Bonnaud, unpub. data). Ac¬ 
cordingly, when comparing the 3’ end of 16S rDNA gene 
sequences, in terms of % substitution, the difference 
between two distant populations of S. officinalis 
(Mediterranean and English Channel) is not significant 
(Table 3). Therefore, this gene portion is not appropriate to 
provide evidence of differences among S. officinalis popu¬ 
lations. 

At higher taxonomic levels, the status of many 
cephalopod groups is still uncertain and needs to be recon¬ 
sidered with the help of phylogenetic reconstruction. 
Phylogenetic hypotheses to estimate divergence time are 
usually given through paleontology (Doyle et al ., 1994). 
True cuttlebones are only known from the early Tertiary; 
Teuthoidea and Sepioidea are thus supposed to have 
diverged in the Cenozoic. 

The use of mtDNA 16S for phylogenetic recon¬ 
struction was shown to be appropriate to test Tertiary diver¬ 
gences in insects and vertebrates (Simon et al, 1990; Hillis 
and Dixon, 1991), but such is not the case for cephalopods 
(Bonnaud et al ., 1994). The phylogeny of cephalopods is 
unresolved at the suprafamilial level because of excessive 
nucleotide divergence (saturation), perhaps due to earlier 
emergence than Cenozoic, or to unequal evolutionary rates 
among taxa. 

To compare results derived from protein elec- 



Octopus vulgaris 
Eledone moschata 
Loligo vulgaris 
Sepia officinalis 
S. orbignyana 
///ex argentinus 
Nautilus macromphalus 
Keyhole limpet 


Fig. 3. Protein (pg/ml) Immunoreactivity of hemocyanin of various cephalopod species against keyhole limpet hemocyanin antiserum, tested by the ELISA 
technique. 






















BOUCHER-RODONI AND BONNAUD: BIOCHEMICAL AND MOLECULAR PHYLOGENIES 83 


Table 3. Nucleotide divergence of 3’ end of 16S rDNA gene sequences at 
various taxonomic levels: (a) populations, (b) intrageneric, (c) inter¬ 
generic. 

Species Divergence 

(% substitution + gaps) 


(a) 

Sepia officinalis (Roscoff) - 



S. officinalis (Banyuls) 

ca. 1 

(b) 

S. officinalis - S. orbignyana 

13.8 


S. orbignyana - S. elegans 

8.1 


S. elliptic a - S. pharaonis 

10.3 


S. pharaonis - S. smithi 

8.8 


Loligo forbesi - L. vulgaris 

9.4 


Nautilus macromphalus - N. pompilius 

5.5 

(c) 

Sepietta sp. - Sepiola atlantica 

5.3 


Sepietta sp. - Rossia macrosoma 

13.0 


Sepiola atlantica - Rossia macrosoma 

12 1 


Sepietta sp. - Sepia officinalis 

17.9 


Sepietta sp. - Loligo vulgaris 

18.4* 


♦Overestimated value, because of an insertion of ca. 20 bases in Loligo. 


trophoresis and sequence data from the 3’ end of 16S 
rDNA, the sequences corresponding to some of the taxa 
appearing in Fig. 1 were aligned and analyzed by 
Neighbor-Joining and PAUP methods. Both approaches 
show that nucleotide analysis clearly separates octopods 
and decapods, but again the monophyly of the order 
Sepioidea, including sepiids and sepiolids is not supported, 
the sepiolids being excluded from the order (Fig. 4). In 
terms of % substitution, sequence comparison shows that 
sepiolids are as distant from Sepia (17.9%) as from Loligo 
(18.4%) (Table 3). Intrageneric divergence ranges from 8 to 
14% for sepiids (the highest value concerning S. officinalis 
and S. orbignyana ), but a direct relationship between taxo¬ 
nomic rank and nucleotide divergence cannot be estab¬ 
lished, as the nucleotide divergence level in other taxa 
could be in another range. Between sepiolid genera, for 
instance (Table 3), the intergeneric divergence is 5.3% 
between Sepiola and Sepietta, two morphologically closely 
related species, but can be as high as 13.0% between 
Sepietta and Rossia. One surprising result concerns the 
sequence comparison of two morphologically distinct 
species of Nautilus which display a low nucleotide diver¬ 
gence value (5.5%). 

DISCUSSION AND CONCLUSION 

The reliability of eye lens proteins as a taxonomic 
tool was shown by some authors (Smith, 1969; Swanborn, 
1971; Brahma and Lancieri, 1979; Tranvouez and Boucher- 
Rodoni, 1990). It is confirmed here that eye lens protein 
electrophoresis analysis serves to group closely related taxa 
together, but its use for inferring phylogenetic relationships 


leads to more questionable results. The use of immunologi¬ 
cal properties of eye lens proteins and of hemocyanin to 
analyze taxonomic relationships is rather disappointing, 
mainly because the distances measured by the ELISA tech¬ 
nique between taxa are not symetrical, /. e. the distance 
between Sepia and Nautilus is not the same as the distance 
between Nautilus and Sepia. The immunological distance 
indicates the degree of similarity, but it is not a character 
that is easy to precisely quantify for phenetic or phyloge¬ 
netic analyses. 

The development of molecular biology has raised 
great hope and excitement about phylogenetic reconstruc¬ 
tions. We now have direct access to the genetic material. 
However the enormity of the available genetic information 
is itself a problem: where is the most appropriate place to 
investigate to answer our questions? Rates of evolution dif¬ 
fer from one gene to the other, and the main problem is to 
find a genetic marker appropriate for the hierarchical level 
we are interested in. If a gene has evolved too rapidly, it 
will be saturated with substitutions and provide a non-sig¬ 
nificant result. If it is too stable, the variability will not be 
informative enough. Our knowledge of genetic structures 
comes primarily from results obtained with vertebrates or 
Drosophila , but their levels of nucleotide divergence are not 
necessarily adequate for analysis of other phylogenetic 
relationships. This cannot be known a priori , and a series 
of prerequisites are necessary before starting a molecular 
analysis. (1) Choice of the gene: its variability must be ade¬ 
quate to the hierarchical level being considered. (2) Choice 
of the species (number and “quality”): this is important 
because some cephalopod genera include over 100 species, 
whereas others are monotypic. (3) Length of the sequence: 
apparently, bootstrap values increase with length of 
sequence, but the reliability depends rather on the number 
of taxa (Lecointre et al. , 1993). So, as with morphology, we 
will have to increase, as much as possible, the number of 
species for each taxon. This is obviously impossible when 
some genera or even families are monospecific. (4) Choice 
of an outgroup: it should comprise more than one species of 
a monophyletic group, close enough to the investigated 
taxon to preclude saturation, but sufficiently different to 
prevent inclusion in the ingroup. 

As far as cephalopods are concerned, nucleotide 
sequence analysis provides a more reliable picture than 
electrophoresis and immunology which are interesting tax¬ 
onomic tools, but are not satisfactory for phylogenetic 
analyses. However, a direct relationship between taxonomic 
rank and % nucleotide divergence cannot be established, as 
the nucleotide divergence level is different in different taxa. 

The congruence between morphological and molec¬ 
ular analyses, two independent sets of data, is very impor¬ 
tant in construction of evolutionary patterns. Morpho¬ 
logical-molecular comparisons however are still rare, and 
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- Eledone cirrhosa 

- Octopus cyanea 

- Rossia macrosoma 

- Sepiola atlantica 

Todarodes sp. 

— Loligo vulgaris 

- Sepioteuthis lessoniana 
- Sepia officinalis 

- Sepia orbignyana 

- Sepia latimanus 

Fig. 4. Phylogeny inferred by Neighbor-Joining distance analysis from mtDNA 16S sequences (Bonnaud et al , 1994). Bootstrap values < 50 are indicated. 


76.5 


100 


50 


66.6 


85 


not often congruent. In cephalopods, one solid and congru¬ 
ent result concerns sepiolids. Their taxonomic and phyloge¬ 
netic position has been a matter of much discussion, being 
classified either as a family of the Sepioidea (Naef, 1912; 
Voss, 1977; Mangold and Portman, 1989) or as a family of 
the Myopsida together with the Sepiidae (Berthold and 
Engeser, 1987). Fioroni (1981, based on embryology) and 
Clarke (1988, based on morphology) proposed to raise the 
sepiolids to ordinal rank. The present electrophoretic, 
immunological, and molecular results confirm that sepi¬ 
olids can be separated from the Sepioidea. 

The analysis of phylogenetic relationships among 
coleoids must take into account a number of difficulties, 
whatever the criteria used, morphological, biochemical, or 
molecular: (1) poor fossil remains of recent taxa; (2) no 
outgroup: we choose octopods as an outgroup for decapods, 
and vice versa because, even if the distance is very impor¬ 
tant, we have no real alternative; (3) many monospecific 
genera; (4) traditional classification not well stabilized; (5) 
little information on rate and modalities of molecular evo¬ 
lution: the evolutionary rate of the different genes is 
inferred mainly from vertebrate results - vertebrates repre¬ 
sent only one episode of the saga of life and it is not always 
possible to transpose to invertebrates genetic postulates 
based on vertebrate results. 

To better understand higher hierarchical levels of 
phylogeny, it appears necessary in many cases to consider 
more than one gene, preferably by associating mitochondri¬ 
al and nuclear information. 
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